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SUMMARY 



Unsteady lift functions for wings of finite aspect 
ratio have "been calculated by approximate methods involv- 
ing corrections of the aerodynamic inertia and of the an- 
gle of the infinite wing. 

The starting lift of the finite wing is found to he 
only slightly less than that of the infinite wing; whereas 
the final lift may be considerably less.. The calculations 
indicate that the distribution of lift near the start is 
similar to the final distribution. 

Both the indicial and the oscillating lift functions 
are given. Approximate operational equivalents of the 
functions have been devised to facilitate the calculation 
of lift under various conditions of motion. 



INTRODUCTION 



The two-dimensional potential theory of airfoils in 
unsteady motion was set forth by Wagner (reference 1) and 
has been extended to problems involving the motion of 
hinged or flexible airfoils by Theodorsen (reference 3) 
and Kiissner (reference 3)., 

It is known that, in the case of steady motion, a 
correction is necessary before the results of the two-di- 
mensional theory can be applied to actual wings of finite 
aspect ratio. No corresponding correction has, however, 
been developed for the conditions of nonuniform motion.* 
The calculations contained herein lead to such a correc- 
tion, which is applicable, within the usual limitations, 
to arbitrary dynamical conditions of motion, 

*A discussion in general terms of such a correction has 
been given in a recent paper by Cicala (reference 4) . 
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THE .JUDICIAL LIFT 
Influence of the Wake 



Owing to the presence .of circulation, the lifting 
wing leaves in its wake a surface of discontinuity whose 
local vortex strength will be determined "by the rate of 
change of circulation taken ho th across the span and along 
the flight path, (See fig. 1.) The formation of such a 
wake and the configuration of the vortices in the wake are 
determined by the assumption that the flow at each instant 
conforms to the Kutta condition. An essential feature of 
the problem is the determination of the influence -of this 
wake on the wing. 

It is important to note that the wake is supposed to 
remain plane and undistorted. The consequence of this as- 
sumption is that the effects' of the wake are additive, per- 
mitting the various flows to be built up by superposition. 

Thus, if the solution is known for the case of a sud- 
den start of the motion - or what amounts- to the same 
thing under 1 the assumptions involved, a sudden change in 
angle of attack - this solution may be used as the element 
in an integral that gives the lift in any variable motion, 
fith this point in mind, attention will first be directed 
to the special case in which the wing starts suddenly from 
rest at t = 0 and maintains the uniform flight velocity 
IT with the normal velocity w thereafter. 

Before the three-dimensional problem is considered, 
if will be helpful to review certain aspects of the two- 
dimension'al theory. Figure 2 shows the elemental flow 
used as a starting po int "by Wagner (reference 1). 

This flow is caused by two vortices, representing an 
element of circulation around the wing and the vortex left 
in the wake when this circulation was generated. The 
streamlines of this flow are eccentric circles and' one 
such circle of unit radius is chosen to represent the wing 
section. The axes are then placed so that this circle has 
its center at the origin. The geometry of the pattern is 
such that, when the wake vortex is at z , ' the wing vortex 
must be at l/z. Such spacing preserves the unit circle as 
a streamline of the flow. 
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Tran'sf prmation of the pattern "by the formula 

H = z + \ ... CD 

flattens the un i ;t circle into a thin-line wing section 
and distorts the originally circular streamlines into oval 
Joukowski figures (fig. 2) 

The transformed pattern thus represents the circula- 
tory flow about a series of symmetrical J.o'ukowski airfoils 
with' an associated count ervort ex in the wake. For conven- 
ience, the ideally thin airfoil is chosen./ 

Each elementary flow of the type shown (fig. 2) con- 
tributes a certain velocity around the trailing edge of 
the airfoil. An instantaneous flow due to an angle of at- 
tack of the airfoil also gives a trailing-edge velocity, 
"but of opposite sense. On this basis,, the problem of cir- 
culation with varying angle of attack may be solved by an 
inverse procedure. Assume some convenient distribution of 
wake vort icity and calculate the trailing-edge velocity at 
each point along the flight path corresponding to the for- 
mation of the proscribed wake. The particular variation 
of angle of attack necessary to cancel this trailing-edge 
flow at each instant (Eutta condition) can then be deter- 
mined. If a number of such curves are found, they may be 
added in various proportions so as to approximate any pre- 
scribed variation of angle of attack; the corresponding 
curves of variation of circulation along the flight path 
are then added in like proportions. 

It was in essentially the manner just described that 
fagner (reference 1) built up the flow pr'odueed by a wing 
during the uniform motion following a sudden unit change 
in angle of attack. The integrated" pressures over the 
airfoil give a lift coefficient that approaches asymptoti- 
cally the known steady value 2tt,' The starting lift is 
found to be exactly one-half this value. The center of 
pressure remains at the quarter-chord point throughout the 
motion. 

In the case of the finite wing, an element of the 
wake will be as depicted in figure 2 but will, in addition, 
contain vortices completing each circuit to the wing 
through the tips,, Hear the start of the motion, these tip 
vortices will be short and their influence on the wing will 
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consequently- "be small, . -Hence , the starting lift of the 
finite wing will be very nearly equal to that of the in- 
finite wing. As the wake lengthens, the correction will 
increase and finally approach the magnitude given by the 
Prandtl theory. 

In the present problem, it is, desired to follow along 
the lines indicated by the Prandtl theory insofar as pos- 
sible , using the existing two-dimensional theory as a basis 
and determining the effect of aspect ratio as a correction, 
fith long tip vortices; and with smooth distributions of 
load, calculations show nearly uniform distributions of 
downwash over the wing. With a short wake, however, the 
calculations show proportionately much greater curvature 
of the induced flow so that the effects near the start 
cannot be approximated by a simple angle-of -at tack correc- 
tion as in the Prandtl theory. 



Lift Hear the Start 

The starting lift of any wi.ng may be* expressed in a' 
very simple manner based on the Kutta condition. It is 
seen that, as a con sequence . of this condition, the portion 
of wake adjacent to the trai ling edge . must move as an im~ 
permeable extension of the wing surface. The flow pro- - 
duc.ed at the first instant is the same as that caused by: 
the wing in process • of growing wider at the rate V. It 
follows that the starting lift -may be thought of as the ■ 
reaction to uniform motion of a wing with increasing mass: 



(2) 



where m 1 is the mass representing the aerodynamic iner- 
tia of the wing. 

. - In order to apply equation (2) to the finite wing, 
it is necessary to. know the inertia factor for such a 
wing as a function of the width. In general, it would 
therefore be necessary to determine the potential flow' 
for normal motion of each wing of given shape. Solutions 
are provided, however, by classical hydro dynamic the'ory 
fo-r. elliptic plates, and it is possible to use these solu- 
tions to represent approximately the initial rate of in- 
crease, of inertia of a wing of conventional shape. 
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The distribution of potential over each chordwi se 
section of an elliptic plate in nor rna 1 motion ha s the 
same form as the corresponding two-dimensional potential, 
At a certain value of the normal velocity, this distri- 
bution is given by a circle with the chord of the section 
as the diameter. Thus 



In two-dimensional flow, the normal velocity is 1 . 00 and 
it is slightly greater for the. finite disk. The factor" 
of increase is the ratio of the semiper ime t er to the span, 
given by the elliptic integral S , i.e. , 



cp = i yrvp 



If the edge of the plate distorts into a slightly 
wider ellipse (fig. 3), the potential differences arising 
will be measured by the differences between the ordinates 
of the original and of the slightly expanded circles, The 
change in the factor S may be neglected. The pressure 
differential at any point is given by the formula 



p » -*2p 



7 22 + ^2 

8x at 



(3) 



where 



cp = 1 ~ x 2 = sin 6 



Bcp 



* x JTT7 



= - cot e 



(4)' 



and, from the geometry of the circle, 
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The pressure across the plat e-i with .normal, velocity w == S 
'and f light velocity IT iW .theref.o : r e f 

p. = PV [2 COt 6 - COt i 1 (6) 
t-o (_ 2 J- 

Intcgration of this pressure oyer any section gives the 
lift coefficient 

o L = % (?) 

with each local center of pressure at the quarter-chord" 
line. 

The start of the elliptic wing is equivalent to a 
uniform lengthening of each chord, so that the true ellip- 
tic outline" is not preserved. It may he shown, however, 
that such a change does conform very nearly to a change 
■into another slightly larger ellipse at all points except 
those near the tips. Furthermore, if it is assumed that 
the wing is distorted in any of a number of ways into a 
slightly different elliptical plan form, it is found that 
the change of virtual inertia is hut little affected "by 
the change in shape and depends primarily on the over-all 
change in size. Each such distortion can be thought of as 
representing a certain di stribution of the starting veloc- 
ity V around the edge of the wing. Equation (7) is exact 
for all these distributions and, since they may be made to 
fall closely on either side of the curve V = constant 
(representing the start of a rigid wing), the equation is 
considered applicable to this case also. 



The Downwash Correction 

A fairly accurate curve of the growth of lift might 
be draTTii by connecting the starting value given by equa- 
tion (7) asymptotically to the known steady value. After 
the wing has progressed several chord lengths, however, 
the effect of the vortex wake can be treated simply as a 
modification of the angle of attack and in this way it 
will be possible to obtain a closer approach' to the true 
form of the curve. 

Figure 4 shows how an elementary loop vortex in the 
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wake of the finite wing' can "be formed by cancelation from 
an element of the wake of the infinite wing, The flow 
produced by- segment s BCD &FH account s for the aspect-ratio 
effect . 

The relation determining the spacing of the vortices 
A and E (fig* 4) is found in the two-dimensional theory. 
(See fig. 2.) At the instant an element of circulation is 
acquired by the wing, the bound vortex and the counter wake 
vortex coincide at the point corresponding to the trailing 
edge of the airfoil. The wake vortex is subsequently car- 
ried downstream with the fluid at- the velocity V while 
the bound vortex moves ahead into the wing as shown. It 
is important to note that the position^of the centroid of 
the wing circulation is unchanged by the transformation' 
from the circular to the flat wing sections. Thus, when 
the wake vortex is at the flattened bound vortex re- 

tains its centroid at l/z. The equivalent length x of 
the tip vortices, in terms of the flight-path distance s, 
become s - 



1 4- s + a/s 3 + 2s - 



1 + s + JT^ + 2 s. 



*= 7s (s + 2) ( 8) 



The tip vortices do no t , of cour se , terminate abruptly 
but merge into the wing. The effective length x gives 
a good approximation to the effect at some distance. Fig- 
ure 5 illustrates the rapid spread of the discontinuity 
into the wing subsequent to its origin at the trailing 
edge. 

Since the vortex wake is supposed to remain plane and 
undi started , the principle of superposition can be applied 
■to the calculation of downwash w i (s). Then the downwash 
due to any distribution of vorticity along the flight path 
and- over the span can be computed by integration of the 
effects of elementary lines of the pattern indicated in 
figure 4. In the simplest aspect of the problem, the sur- 
face of discontinuity may be replaced by a "skeleton' 1 com- 
posed of two tip vortices of finite strength connected * 
across the span by a vortex sheet. The tip vortices are 
located in such a way that their influence at the middle 
of the wing approximates the inf luence of the spanwi se 
grading of discontinuity. The problem will first be treat- 
ed in this manner and will later be extended to include 
the elliptical form of spanwi se loading. 
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A straight down wash flow at the center "being" assumed , 
the calculation for an elementary loop is carried out at 
the chordwise cent r old of the wing vorticity where the 
small influence of segments 33 and & disappears. By 
the application of Biot-Savart 1 s rule, the down wash due to 
segments DO and IE is found to he 




The transformation of the length x to the distance s 
is accomplished "by equation (8). 

With the tip' vortices at the fixed spacing y, the 
downwash due to any variation of vorticity along the flight 
path 7 ( s ) may he found from 

dw^ f? 
Wi(s) - — (a) 7(0) + / — (s - b q ) 7 f (s ) ds 0 (10) 
d7 7 ^7 0 0 0 

J o • 

Circulation and Lift 



The circulation around the wing at later stages of" 
the motion may bo determined from the two-dimensional 
theory by using the effective angle of attack 

a e = a - aj (n) 

where 

w i 

a ± = -~ 

v 

Let l\ (s) be the rise of circulation following a unit 
jump of angle of at tack (indi cial circulation) as given 
by Wagner (reference 1) for the infinite wing. (The sub- 
script i is used hereinafter to denote the effect of a 
unit jump of angle of attack occurring at s == 0.) Then, 
for the finite wing, 



7(b) « r x (s) o e <0) + /" P r (s - s 0 ) ae 1 (s 0 ) ds 0 (12) 

4 



-ST., A * C . A . Teehrti cal ■ No t e , Hp V . 682 



9 



The determination of the. eff^ctiye angle of attack and the 
circulation of the finite wing thus, depend on the simul~ 
taneous solution of integral equations (10) ^.nd (12). 

The formal solution of equations (10) , (11), and (12) 

for the indicial downwash .w* (s) and for the indicia! 

x 

circulati on y x (s) would be expected to he quite diffi- 
cult* It is, however, a fairly easy matter to find satis^ 
factory curves by trial, particularly since it is known 
that 7i(s) coincides with ? x (s) and that w^ (s) =0 
at the beginning (s = 0). Figures 6 and 7 show curves 
determined in this way for two different values of y, 

corresponding to aspect ratios 3 and 6 (y ~ £ for the 

tt £ 

elliptical wing). 

The lift at the later stages of the motion is found 
by combining the effective angle of attack a e : (s) with 
the two-dimensional indicial-lif t function c^ (g) of 

Wagner. The total indicial lift for the three-dimensional 
wingy , neglecting the small increment due . to acceleration 
of the downwash flow, is therefore 

C Ll (s) = c, (s) -h / c l-(s - s 0 ) ae x * (s 0 )ds 0 (13) 

^o 

Plots of Ct ji determined by graphical integration are 

given in figure 8, . 

The lift of the finite wing approaches a steady value 
much more rapidly than indicated by Wagner's curve. It is 
observed that the downwash due to the tip vortices in- 
creases so slowly that an effective reduction of the angle 
of attack does not occur until the lift has risen nearly 
to its asymptote. 



Lift with Elliptical Span Loading 

The foregoing calculations of 0Ce(s) and T(s) were 
applied to wings with finite tip vortices. Any preassigned 
variation of the loading across the span may be built up, 
however, by the addition of elementary vortex loops. 
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The expression for- the downwash at the center of the 
wing may be integrated in the case of elliptic loading. 
Let 

7 = 7 0 sin 6 (14) 



17 = 7 Q cos e d'6 



(15) 



and let y = — cos 6 
2 



y 2 + X s 



-[(*)' -It 



1 - k 2 sin 2 6 



] 



where 



k = 




(16) 



Then the downwash due to a series of loops of the form 
CEF (fig, 4) is 



O TT/2 



W = y 

2tt C 



w = — 7 

2tt 



I cos 2 e \ 



d e 



(17) 



•k 2 sin 8 G 



° (I)' 



rk - . K (k) + i 



(18) 



where K(k) and l(k) are the complete elliptic inte- 
grals. (See peirce ! s table.) 

The induced downwash per unit circulation around the 
wing., as defined under The Downwash Correction, is 



dw* 1 

— (*) = — 

d7 2tt 



xk 



— X(k) + - 

(I) 2 



K(k) (k - i) + lE(k) - 1 1(19) 
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It is necessary, as before, to find the equivalent 
length x from the distance s V (See equation (8) . ) 

The indicial-lif t functions (fig. 8) were determined 
by graphical integrations, as in the preceding case. 



OPERATIONAL F OEMULAS FOB THE LIFT 
Indicial Lift Functions 



The lift of a wing under various dynamical condi ti ons 
may "be conveni ent ly found by the operat i onal method de- 
scribed in reference 5, In order to facilitate the use of 
the lift functions in such problems, approximate formulas, 
together with their operational equivalents, have been de- 
vised; 

C L (s) = C 0 + C i e r ^ s + C a e r 2 s (20) 
Values of the constants are listed in the following table-: 



TABLE I 



A 


Co 








r 2 






-0.330tt 


-0.0455 


-0. 670tt 


-0.300 


6 


4.71 


-1.740 


- . 324 


0 


0 


3 


3.77 


-1.07 


-.490 


0 


0 



Points shown in figure 8 illustrate the degree of ex 
actne s s of these formulas . 



The operational equivalents of the preceding expres- 
sions are readily formed from the relation 

e rs = — — l(s) (21) 
D - r 

whence ; 

C L (s) = C L (D) l( s ) = (C 0 + C x — + C a \X(s) (23) 

1 ^ D - r x D - r 2 / 
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Lift 1 Functions for an Oscillating Airfoil 
Poll owing the general formula 

C L (s) = C Li (D) a(s) (23) 
the lift in sinusoidal motion, where 

■ a - e lns (24) 

is given by 



C T (s) = C L (D)e ins * = C Tj (D) .,, P > l( s ) (25) 



X) - in 

where the instantaneous lift due to acceleration is omit- 
ted. On expansion of thi s operator it is found that , 
with the exception of transient terms, 

C L (s) = C L (in) e ins (26) 

The real and the imaginary " part s of Cj^ ( in) corr e- 
spond to the functions \F(n) and G-(n) developed for 
the two-dimensional case by Theodorsen (reference ; 2) : 



C Ll (in) = 2rr [F(n) + iG(n)] (27) 



where 



2 2 

n n 



2ttP = C 0 -h C x - — + C 



r x 2 + n s r 2 s + n 2 



2rr& = -C a ^ C 2 &g 



r 1 3 + n 2 , r 2 2 + n 2 

Figure 9 shows these functions as calculated from the val- 
ues listed in table I. 



Langley Memorial Aeronautical Laboratory , 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , December 14, 1938. 
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Figure 6.- Indie ial-downwash curves. 



18 20 22 24 

















1 1 i i i i i i 
Finite tip vortices 


































- — --Elliptic loading 
























































A 


sped ral 


Ao 


































































Wa 


-gn< 


?r's 


c 


jrv 


























































































































6" 




































































































"3- 









































































































































































































































































































































































s,half chords 
Figure 7,- Indicial- circulation curves. 
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